To Appear in ApJL 

Preprint typeset using 1^1^^ style emulatcapj v. 8/13/10 



O 
(N 

> 

o 
O 

00 



I 

o 



> 
in 
o 

m 
d 

o 



A POPULATION OF VERY-HOT SUPER-EARTHS IN MULTIPLE-PLANET SYSTEMS SHOULD BE 

UNCOVERED BY KEPLER 

Kevin C. Schlaufman"''^, D. N. C. Lin^'^, and S. Ida^ 

Received 2010 July 22; accepted 2010 October 12 

ABSTRACT 

We simulate a Kepler-\\k.e observation of a theoretical exoplanet population and we show that the 
observed orbital period distribution of the Kepler giant planet candidates is best matched by an 
average stellar specific dissipation function Q'^ in the interval 10^ < Q'^ < 10''. In that situation, 
the few super-Earths that are driven to orbital periods P < 1 day by dynamical interactions in 
multiple-planet systems will survive tidal disruption for a significant fraction of the main-sequence 
lifetimes of their stellar hosts. Consequently, though these very-hot super-Earths are not characteristic 
of the overall super-Earth population, their substantial transit probability implies that they should 
be significant contributors to the full super-Earth population uncovered by Kepler. As a result, the 
CoRoT-7 system may be the first representative of a population of very-hot super-Earths that we 
suggest should be found in multiple-planet systems preferentially orbiting the least-dissipative stellar 
hosts in the Kepler sample. 

Subject headings: planetary systems — planets and satellites: formation — planets and satellites: 
individual: CoRoT-7 — planet-star interactions 



1. INTRODUCTION 



NASA's Kepler mission (iKoch et all 120101: 
Jenkins et a l.ll 2010al: [Caldwell et al.ll2010t iBrvson et al l 
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Batalha e t al. | l2010al |bl: (Haas et al J 120101 : 

is currently searching for transiting 
planets in photometric observations of ~ 156, 000 
solar-type stars in a 115 deg^ field of view in the 
constellation Cygnus toward the Orion arm of our 
Galaxy. Kepler photometric observations will achieve 
80 parts per million precision over 80%-90% of the 
telescope field of view, sufficient to identify the transit 
signal from an Earth-Sun equivalent system. By the 
end of its four-year mission, the combination of Kepler^s 
unprecedented precision, sample size, and homogeneous 
selection will produce a near ideal set of exoplanet 
detections to compare with theoretical models of planet 
formation. Though the preliminary exoplanet candidate 
list announced in Borucki et al. (2010) is incomplete, 
it is still a valuable constrain on theoretical models of 
close-in planet formation. 

In the core accretion model of close-in pla net forma- 
tion fe.g.. lPollack"etaI1ll996l:rida fc Lin|[200l . the cores 
of giant planets form near the ice line of their parent pro- 
toplanetary disk. The cores grow to an isolation mass 
and accrete gas until they are massive enough to open- 
up a gap in their parent disk. The newly-formed giant 
planets then Type II migrate into the close proximity of 
their host star and stop inside the magnetospheric trun- 
cation radius o f their parent protoplanetary disk (e.g., 
iLin et al.l[T99l . Lower -mass planets in the close prox- 
imity of their host stars result from the inward Type I 
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migration (|Wardlll997[ ) of planetary embryos that did not 
reach the mass necessary to initiate runaway accretion of 
gas and become giant planets. These lower-mass close-in 
planets may be hot Neptunes or icy super-Earths that 
formed near the ice line of their parent protoplanetary 
disks and then migrated to just oi itside the disk mag - 
netospheric truncation radius (e.g.. iMasset et alll2006[ ). 
Alternatively, they may be rocky Earth and super-Earth 
mass planets that formed interior to the ice line and 
then migrated either singly or in resonant planet con- 
voys into the close proximity of their host star, eventu- 
ally stopping just outside the disk magnetospheric trun - 
cation radius (e.g.. llda fc Lh]l2010l : lOgihara et al.ll2010( ). 
Theoretical exoplanet population synthesis (EPS) mod- 
els (Ida & Lin 2008, 2010; Mordasini et al. 2009a^b) in 
the core accretion and migration paradigm have repro- 
duced many trends in the observed distribution of exo- 
planet properties, and comparisons between observations 
and models are especially useful to constrain u ncertain 
model parameters (e.g., Schlauf man et al]l2009l ). 

Tidal interactions between close-in planets and their 
stellar hosts contribute significantly to the dynamical 
evolution of close-in planet systems after th e dissipation 
of their parent protoplanetar y disks (e.g., iRasio et"al] 
[T99llDobbs-Dixon et al.1[200l . Tidal effects tend to cir- 
cularize and shrink orbits to the point that close-in plan- 
ets suffer sub stantial rnass loss and eventually tidal dis- 
ruption (e.g.. [Gu et al.l l200l iJackson et aHl2008al I200I 



ILi et al.ll201C ). The efficiency an d timescale for these 
processes are both uncertain (e.g.. ILevrard et al1l2009D 
and highly depende nt on th e internal structure of the 
close-in pl anet (e.g.. IQgilvie fc Lin 2004) an d its stellar 
host (e.g., lOgilvie fc LinTEoOTi : iBarker fc Oglivig I200I 
l2010f ). Modeling of the properties of the close-in planet 
candidates identified by Kepler will help resolve some of 
these uncertainties. 

In this lette r, we combine an exten ded version of the 
EPS models of llda fc LinI ([Moh - llLTol hereafter-with the 
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mass-radius relations of iFortnev et al.l (|2007f) to extract 
the expected period-radius distribution of the extended 
EPS. We use Monte Carlo simulations in concert with a 
detailed model of Kepler sensitivity and a simple model 
of tidal evolution to determine the exoplanet period- 
radius distribution that Kepler would observe in the ex- 
tended EPS as a function of the strength of tidal evo- 
lution. We then compare the result with the observed 
period distribution of Keple r transiting gian t planet can- 
didates announced in iBorucki et al.l (|2010 r) to constrain 
the strength of tidal evolution in the Kepler candidate 
planet sample. We describe our Monte Carlo si mula- 
tions and Kepler-like observations of the extended lILld 
EPS in Section 2. We discuss the results and implica- 
tions of our simulations in Section 3, and we summarize 
our findings in Section 4. 

2. ANALYSIS 

Transit observations directly measure exoplanet radii 
and orbital periods, while radial velocity observations are 
necessary to confirm individual transit detections and 
to determine exoplanet masses. Unfortunately, most ex- 
oplanet candidates identified by Kepler will orbit host 
stars with faint apparent magnitudes. As a result, radial 
velocity follow-up will be telescope-time intensive for the 
more massive Kepler exoplanet candidates and possibly 
too imprecise at faint magnitudes to confirm lower-mass 
candidates. Consequently, it becomes important to care- 
fully consider what insight can be gained from the ex- 
oplanet period distribution derived purely from Kepler 
transit photometry. 

To that end, we use a Monte Carlo simulat ion to sys- 
tematically observe an extended version of the|lLi3EPS 
to as much as possible match Kepler observations of 
the Milky Way's exoplanet population. In this case, 
we have extended the ILIO EPS such that the distri- 
bution of host stellar mass in the 10'^ systems in the EPS 
matches the distribution of host stellar mass expected in 
Kepler observations. That is, we matc h the distribution 
of host stellar mass in the extended IILIOI EPS to the 
mass distribution expe cted in a Kepler-like observation 
of the iReid fc HawlevI (,2005) solar neighborhood lumi- 
nosity function, given Kepler^s sensitivity as a function 
of host spectral type and apparent magnitude. 

The first step in our calculation is to appro ximate the 
radius of each planet in the extended lILld EPS. The 
IlLl O models include the asymptotic semimajor axis, rock 
mass, ice mass, and gas mass of every planet. From 
th ose quantities, we use the mass-radius relations given 
in IFortnev et al.l (|2007[ ) to assign radii to giant planets 
as a function of semimajor axis, solid mass, and gas mass 
and to terrestrial planets as a function of solid mass and 
composition. 

We then do 100 iterations of a Monte Carlo simulation 
in which we assign each planetary system in the extended 
ILIO EPS a random host stellar mass, age, and appar- 
ent magnitude, as well as a random orbital inclination. 
We assign each eccentric planet a random argument of 
periastron. We first randomly assign spectral types and 
apparent magnitudes to the host stars of the extended 
exoplanet population. We use the sola r neighborhood 
lumin osity function given in Table 8.3 of IReid fc HawlevI 
(|2005f ) to determine the number density of FGK stars as 
a function of spectral type and apparent magnitude; we 



use this information to determine the probability that 
a randomly selected star from the Kepler field is of a 
given spectral type and apparent magnitude. The spec- 
tral type of the assigned host star determines the host 
stellar mass and radius. 

To account for the effect of tid al evolution on the ob- 
servability of the extended lILld exoplanet population, 
we exploit the fact that the two-Gyr moving-average 
smoothed star formation history in the Milky Way has 
been more or less constant over the last ~ 10 Gyr (e.g., 
iRocha-Pinto ct al. 2000a, b). As a result, the unknown 
age of a randomly selected star from a magnitude-limited 
survey of Milky Way disk stars with main-sequence life- 
time T* short relative to that ~ 10 Gyr interval should be 
distributed more or less uniformly between zero and its 
main-sequence lifetime. Magnitude-limited transit sur- 
veys like Kepler are biased toward stars with Af* > 1 M© , 
and those stars have main-sequence lifetimes r* < 10 
Gyr. Consequently, the unknown system age Tgys of a 
typical candidate exoplanet system identified by Kepler 
in transit with host stellar mass M* should be well- 
approximated as uniformly distributed in the interval 
< Tgys < T*(M*). For that reason, in the Monte Carlo 
we randomly sample the age of each exoplanet system 
from a uniform distribution between zero and the main- 
sequence lifetime of its stellar h ost. For Sun- like stars 
L*(M*) = Lq (M*/Mq)3-^ fe.g.. lPoppeilfl98l . and the 
total amount of hydrogen available for fusion is propor- 
tional to Ma. . The main-sequence lifetime is roughly then 
r,(M*) = tq M,IL,{M,) = 10 (M*/Mo)-2-5 Gyr. We 
consider an exoplanet observable if its randomly deter- 
mined age Tgys is less than the timescale for the tidal evo- 
lution of its orbit to move it within 1 i?* of its stellar host 
Tdis at which point we assume that the planet is tidally 
disru pted and no longer observable (e.g.JSand quist et al.l 
Il998f ). In other words, a planet is only observable if 
Tsvs < Tdis- We define the time to disruption as (e.g., 
Ilbgui fc Burro^l2009f) 
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where ag is the initial semimajor axis of the planet before 
tidal evolution, Q'^ is the specific dissipation function of 
the host star, G is Newton's gravitational constant, Mp 
is the mass of the planet, and R,, is the radius o f the h ost 
star. The expression given in Ilbgui fc Burrow^ (|2009[ ) is 
strictly valid only for circular orbits, so we assume that 
the timescale for eccentricity damping is short relative to 

Tdis- 

As a result of the assumptions we make in our treat- 
ment of tidal evolution, we assume that the all planets 
from the extended JLIQ population with orbital period 
P < 10 days are on circular orbits. We use the eccen- 
tricity distribution generated from the EPS models for 
longer-period systems. We sample the orbital inclina- 
tion and argument of periastron from the standard ran- 
dom distributions of th ose quantities. We ado pt the an- 
alyt ic formulae given in iSeager fc Mallen-Orne las (2003i) 
and lFord et al.l ()2008f ) to determine the transit depth and 
transit duration of every planet found to transit. Finally, 
we use a detailed model of Kepler sensitivity (D. Koch 
et al., private communication) that gives the threshold 
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detectable exoplanet radius as a function of host spectral 
type, host apparent magnitude, transit depth, transit du- 
ration, and number of transits over a given duration of 
observation to determine which of the simulated transit- 
ing exoplanets would be detectable in the first 43 days 
of Kepler data. 

We save the results of this iteration, and repeat our 
Monte Carlo simulation 100 times to average over all 
of the given distributions to generate a prediction for 
the characteristics of the planetary systems Kepler would 
likely ide ntify in an exoplanet population matching the 
extended HLia EPS. We carry out this Monte Carlo simu- 
lation for four fiducial values of Q'^: 10^ (strong tidal evo- 
lution), 10^ (moderate tidal evolution), 10'' (weak tidal 
evolution), and oo (no tidal evolution). We present the 
results of our simulations for each assumed value of Q'^ 
in Figures [T] to m 

3. DISCUSSION 

The exoplanet period-radius distribution t hat re sults 
from a JTepZer-like observation of the extended IILIOI EPS 
is a strong function of Q'^. In the most dissipative case 
Q'^ = lO'"', hot Jupiters are quickly disrupted and the pe- 
riod distribution of giant planets is fiat with a peak at 
P ~ 10 days. At the same time, the handful of planets 
at very short period P < 1 day produced by dynamical 
interactions between planets in multiple-planet systems 
after the dissipation of their parent protoplanetary disk 
are tidally disrupted and unobservable in a Kepler-like 
observation. In the less dissipative cases Q'^ — 10^ and 
Q'^ = lO'^, hot Jupiters survive for a significant fraction 
of the main-sequence lifetime of their host star and the 
period distribution of giant planets peaks in the range 
3 < P < 5 days in agreement with the giant planet pe- 
riod distribution in Borucki ct al. (2010). In those cases, 
the observed period distribution of super-Earths would 
have a long tail toward periods shorter than one day. In 
the case with no tidal dissipation Q'^ = oo, hot Jupiters 
survive for the entire main-sequence lifetime of their stel- 
lar hosts and the period distribution of giant planets is a 
monotonically decreasing function of period. As a result, 
a Kepler-like observation would be domina ted by planets 
with periods less than one day. Unlike the IILIOI models, 
our calculations suggest that in the case of weak or no 
tidal evolution super-Earths should be observed to have 
a shorter average period than giant planets. 

In Figure [5] we compare the giant planet period distri- 
bution that results from our simulated Kepler-like obser- 
vation of the extended iILlQ EPS w ith the Kepler planet- 
candidate period distribution from lBorucki et al]()2010[ ). 
Massive close-in planets are most sensitive to tidal evolu- 
tion, so we only compare t he per iod distribution of giant 
planets from the extended IILIOI EPS more massive than 
50 with those Kepler planet candidates with radii 
greater tha n 0.5 Jupiter radii. T he Kepler planet candi- 
dates from lBorucki et all (|2010D are not necessarily rep- 
resentative of the properties of the still-embargoed full 
sample of Kepler planet candidates. At the same time, 
the announced candidates are very likely those candi- 
dates unsuitable for radial velocity follow-up. For the 
massive planet candidates at least, that only indicates 
they orbit apparently faint host stars. Properties like 
the average metallicity of a stellar population are not a 
strong function of Galactocentric radius, so there is no 



reason to believe that the subsample of planets that or- 
bit apparently faint stars at distances d ~ 500 pc from 
the Sun is systematically different than the subsample 
of planets that orbit apparently bright stars at distances 
d ^ 200 pc from the Sun. For those reasons, we ar- 
gue that a comparison of the period distribution of the 
announced Kepler exoplanet candidates with the period 
distribution expect ed from a Kepler-like observation of 
the extended lILld EPS under different assumptions for 
the strength of tidal evolution is meaningful. 

The period distribution of the Kepler planet candidates 
is inconsistent with both strong tidal evolution and no 
tidal evolution given the extended iILlQ EPS. The ob- 
served distribution is consistent with the two interme- 
diate values of Q'^, = 10^ and Q'^, = lO"^, with a better 
match provided by Q'^, = 10''. This measurement is also 
consistent with the typically quote d Q'^ ~ 10® from the 
literature (e.g.. lOgilvie fc Linll2007l ). 

The average specific dissipation function Q'^ ~ 10'' 
suggested by our analysis implies that though very-hot 
super-Earths are not common, their substantial tran- 
sit probability ensures that they will be readily observ- 
able by Kepler. We expect that about 10% of de- 
tected planets with P < 10 days and Alp < 10 
should be very-hot super-Earths with P < 1 day. Very- 
hot super-Earths are produced by dynamical interac- 
tions in multiple-planet systems, and we suggest that 
very-hot super-Earths systems identified in transit sur- 
veys should frequently have observable companions (e.g., 
iMardling fc"Linll2004i) . Indeed , the very-hot super-Earth 
CoR,oT-7b (ILeger et al.l 120091) is in a confirmed multi- 
ple system (Queloz et al. wit h the rece ntly sug- 
gested possibility of a third planet ( Hatzes et al. 2010). 
Very-hot super-Earths should also preferentially be found 
around Sun-like stars with M* > 1.25 Mq, as these stars 
are likely less dissipativ e than stars with w 1 Mq 
([Barker fc Ogilviel 120091 ) . At the same time, there are 
other explanations for the properties of the CoRoT-7 
syste m that match the observations fe.g.. iJackson et all 
l2OT0h . 

There are many limitations to our approach. We use 
only a single value of Q'^ to model the strength of tidal 
evolution regardless of host star age or spectral type, 
when i s there is e vidence fr om both observation (e.g., 
iSchlaufman 2010: Winn et al. 2010) and theoretical mod- 
els fe.g.- lBarker fc Qgilv ic 2009, 2010) that Qi is a func- 
tion of host star age and stellar structure, as well as 
exoplanet orbital period. This issue will be resolved in 
the future when the full list of Kepler planet candidates is 
announced, as the increased host star statistics will allow 
us to do similar calculations in bins of host star effective 
temperature and thereby constrain Q'^ as a function of 
Tcff and therefore stellar structure. We do not simul- 
taneously evolve the orbit an d radius of the exoplanet 
population (e.g.. [Jackson ~ al. 2008lJ llbgui fc Burrowsl 
12001 IMiller et al.ll20 09V and the migration stopping con- 
diti ons dis cussed in Section 2.5 of ILIO are also uncertain. 

In llLld , planets that were massive enough to open-up a 
gap in their parent protoplanetary disk were stopped just 
inside the disk magnetospheric truncation radius (taken 
to be at P = 2 days), while planets that were not mas- 
sive enough to open-up a gap were stopped just outside 
of the magnetospheric truncation radius (taken to be at 
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P = 3 days). Though these stopping conditions are de- 
pendent on the unce rtain prop erties of the parent proto- 
planetary disk (e.g iLin et al.|[T99l iMasset et al.ll2006l: 
lOgihara et al.l[2010l ). they are self-consistent in that they 
describe where an exoplanet of a given mass will stop if 
the stellar magnetic field and disk mass-accretion rate 
have a given value. 

4. CONCLUSION 

We coupled a detailed model of Kepler sensitivity with 
the results of the extended lILlQ EPS, a simple model 
of tidal evolution, and published mass-radius relations 
for both terrestrial and giant planets to determine the 
perio d-radius distribution expected under the extended 
|IL13EPS. We found that the period distribution of close- 
in giant planets is a strong function of the strength of 
tidal evolution parametrized as Q'^, and we showed that 
the period distribu tion of th e Kep ler planet candidates 
announced by Bor ucki et al.l ()201Clf ) is best matched by 
the extended lILia EPS with 10^ < < 10^. In that 
case, there exists a population of very-hot super-Earths 
with periods less than one day that results from dy- 
namical interactions in multiple-planet systems. Though 
these very-hot super-Earths are rare, their relatively high 
transit probability and long time-scale for tidal disrup- 
tion when 10® < Q'^ < 10^ indicate that they will be ob- 
servable by Kepler. The predicted very-hot super-Earths 
are analogous to the CoRoT-7 system, and we suggest 
that Kepler should find many such very- hot super-Earths 
in multiple-planet systems preferentially around stellar 
hosts with convective cores and radiative envelopes. 
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Figure 1. Simulated Kepler observation of the extended IILICI EPS models assuming the mass-radius relation from lFortnev et al.l lj20071 ) 
and strong tidal evolution (i.e., Q'^ = 10^). The blue shading shows the underlying initial extended EPS before tidal evolution, while the 
blaek contours show the distribution of exoplanet properties including tidal evolution that would have been observable by Kepler in its first 
43 days of science operation. We plot known exoplanets as red squares and confirmed Kepler planets as yellow squares. Left: Period-radius 
plane with marginal distributions. In this case, the orbits of giant planets quickly degrade due to strong tidal evolution. Right: Semimajor 
axis— mass plane with marginal distributions. 
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Figure 2. Same as Figure [Tl but assuming moderate tidal evolution (i.e., Q', = 10^). In this case, the moderate tidal evolution allows 
some giant planets to survive. 



Very-Hot Super-Earths 



7 




Figure 3. Same as Figure [T] but assuming weak tidal evolution (i.e., Q', = 10^). In this case, the weak tidal evolution allows many 
giant planets to survive. Interestingly, the observed period distribution of super-Earths extends to P < 1 day, as super-Earths scattered to 
short-period orbits through dynamical interactions in multiple-planet systems would persist for a significant fraction of the main-sequence 
lifetime of their host stars. These very-hot super-Earths should occur in multiplc-planot systems, and should preferentially orbit the 
least-dissipative host stars in the Kepler sample (i.e., those stars with Af, > 1.25 Mq, convective cores, and radiative envelopes). 
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Figure 4. Same as Figure[T]but assuming no tidal evolution (i.e., Q'^^ = oo). In this case, the lack of tidal evolution allows all very-short 
period systems produced by dynamical interactions to survive. 
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Figure 5. Histograms describing the period distribution of short-period giant planet s. The black c urve represents the period distribution 
of short-period giant planets (e.g., planets with Rp > 0.5 Rj and P < 10 days) frorn Borucki_et al. (2010). The gray rectangles represent 
the one-sigma range in the period distribution of short-period giant planets (e.g., planets with Alp > 50 Mq and P < 10 days) expected 
from a Kepler-like observation of the extended ILIO EPS. Each panel corresponds to different assumptions for the strength of tidal evolution 
in the simulated Kepler observation of the extended ILIQ EPS. Given the extended EPS, the first 43 days of Kepler science data already 
rule-out both strong tidal evolution and no tidal evolution; the current data favors lO" < Q', ?■ 10^. 



